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Modeling of the temperature-dependent spectral response of Ini_ x Ga x Sb infrared 
photodetectors. 


Abstract 

A model of the spectral responsivity of Ini_ x Ga x Sb p-n junction infrared photodetectors has 
been developed. This model is based on calculations of the photogenerated and diffusion currents in 
the device. Expressions for the carrier mobilities, absorption coefficient and normal-incidence 
reflectivity as a function of temperature were derived from extensions made to Adachi and Caughey- 
Thomas models. Contributions from the Auger recombination mechanism, which increase with a rise 
in temperature, have also been considered. The responsivity was evaluated for different doping 
levels, diffusion depths, operating temperatures, and photon energies. Parameters calculated from the 
model were compared with available experimental data, and good agreement was obtained. These 
theoretical calculations help to better understand the electro-optical behavior of Ini_ x Ga x Sb 
photodetectors, and can be utilized for performance enhancement through optimization of the device 
structure. 
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1. INTRODUCTION 


The spectral response is an important characteristic for defining the performance of a 
photodetector. This characteristic specifies the responsivity, defined as the generated photocurrent 
per unit optical power, variation with wavelength. The photodetector responsivity is primarily 
dependent on the device material, structure and the operating conditions in tenns of bias voltage, 
temperature, and wavelength of the incident radiation. It thus becomes crucial to have a thorough 
understanding of the effects of these parameters in order to design and fabricate an optimal 
photodetector. A theoretical model for the spectral response is particularly useful to accomplish such 
task. 

Antimonide based ternary III-V compounds are good materials for photodetectors operating 
in the infrared range. These devices have several applications such as atmospheric remote sensing, 
hazardous gas detection, and optical fiber communication 1 . InGaSb ternary alloy systems show a 
promising performance in the 1.7 to 5 -pm wavelength range. By varying the indium composition, 
the detection wavelength could be tuned to optimize the quantum efficiency. Epitaxial growth of 
InGaSb on different binary substrates has been reported using different techniques". Besides the 
complexity of these techniques, tuning to a specific wavelength usually involves performance 
deterioration mainly due to lattice mismatch problems. Lattice mismatch increases dark current and 
noise, limiting both dynamic range and sensitivity of a detector. The availability of bulk ternary 
substrates significantly simplifies the fabrication process by using simpler and lower cost 
techniques. 3 

In this paper, a model is developed for the responsivity of Ini_ x Ga x Sb p-n photodetectors. The 
model includes the temperature effects on the various parameters and their influence on the spectral 
response of the device. Also the significance of using a ternary substrate will be presented. The 
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model is based on studying the excess minority carrier generation, recombination and diffusion by 
applying the continuity equations. 

The device responsivity, and therefore quantum efficiency, is affected by the absorption 
coefficient, reflection coefficient and minority carriers recombination lifetimes. The absorption and 
reflection coefficients are calculated using Adachi’s model, 4 " 8 modified to include temperature 
effects. The mobility is calculated using an extended Caughey-Thomas model, 9 ' 11 which accounts for 
both temperature and carrier concentration. The main mechanisms acting on the minority carriers 
lifetimes are the Shockley-Read-Hall (SRH) recombination, 12 " 14 and Auger recombination 15 " 19 
besides the surface recombination effects. The Auger processes are strongly dependent on 
temperature at high concentrations and have a limiting effect on the device responsivity. Validation 
of the model is performed through comparison with experimental characterization of actual 

Ino. 2 Gao. 8 Sb photodiodes. These photodiodes were manufactured by simple diffusion technique, 

20 

using an InGaSb ternary substrate." 

The paper is organized as follows; in section 2, the methodology is introduced for calculating 
the temperature dependency of the energy bands, optical dispersion relations, and the effects of 
mobility and recombination lifetime on diffusion length, all of which are required for determining 
the photocurrent densities and spectral responsivity using the continuity equations. In section 3, 
results obtained from the model are presented and discussed, giving emphasis to the variation of the 
responsivity with temperature and the influence of using the ternary substrate. Finally, conclusions 
are given in section 4. 
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2. METHODOLOGY 


The modeling of the thermal characteristics of the spectral response is based on calculating 
the device parameters considering their relationship to the energy band structure of the material. The 
schematic for the modeled photodiode structure is shown in Fig. 1. These photodiodes were 
manufactured by Zn diffusion to form the p-layer on an n-type InGaSb ternary substrate. 20 The 
resulted photodetector has a p-n junction structure illuminated from the top p-layer. Considering a 
non-degenerated semiconductor with full ionization, the model assumes a steady-state solution, with 
low electric field not to saturate the carriers velocity. 

2,1 Lattice parameters and energy bands 

Ini_ x Ga z Sb is a ternary alloy whose parameters and characteristics change according to its 

composition, %. The lattice parameters for Ini_ x Ga x Sb are obtained from the linear interpolation of 

2 1 

the transition parameters of the binary materials InSb and GaSb using the relation" 

Qalloy = X'^GaSb + 0 — X)^lnSb _ XO ~ X )C B 5 (0 

where B and Q represent the binary and ternary parameters, respectively, and C« is the bowing 

8 21 

parameter that accounts for contributions arising from lattice disorders. ’ 

Applying Eq. 1 , the electron effective mass corresponding to the conduction band, m e , heavy- 
hole band, nn„ and spin split-off band, m so , as well as the relative permittivity, s s , for Ini_ x Ga x Sb, are 
given by m e = (o. 0 1 5 + 0. 0 1 % + 0. 025 j 2 ) • m o , m h = (0. 43 - 0. 03 x) • m 0 , m so = (0. 1 9 + 0. 05 %) • m 0 , 

and s s = (l 6.8 - 1 . 1/)- , respectively, where m„ is the electron rest mass, and s 0 is the absolute 

• • • 22 
permittivity. 
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The energy bandgap, E, variation with temperature, T, due to the temperature-dependent 
dilation of the lattice and electron-lattice interaction, is approximated by the semi-empirical Varshni 
equation 23 

E(T) = E(0) — , (2) 

(T + J3) 

where E(0) is the energy bandgap value at absolute zero, and 8 and f:i are fitting constants. Table I 
lists GaSb and InSb fitting constants for the different critical-point energies. Parameters for E a , Ej, 
Ej+A i, and E 2 were taken from reported experimental data," ' while the E a +A 0 temperature 
dependency for InSb was assumed to be parallel to E 0 ’ . Generally, the energy bands present an 
almost linear temperature dependence, with gradual increase at lower temperatures. Once the 
temperature dependent critical-point energies are defined for the binary materials, eq. (1) is applied 
to find the corresponding ternary material energies and its temperature dependence, with C« = 0.2 eV 
only for fs 0 . : 8,21 

2.2 Optical dispersion parameters 

The optical properties of a semiconductor are described by the complex dielectric function 
£(X,T ) = s x (X,T ) + i ■ £ 2 (X,T ) , at all photon energies he/ A , , where h is Plank’s constant, c is the 
speed of light and X is the radiation wavelength. The real and imaginary parts of the dielectric 
function, Si(X,T) and £ 2 ( 8 , T), respectively, are analytically related using the Kramers-Kronig 
relation. The imaginary part is calculated from the joint density of states function at various critical 

A 8 

point energies in the Brillouin zone and indirect band-gap transitions, as shown by S. Adachi. ' 
Thus, the complex dielectric function is strongly related to the energy-band structure of the medium 
and its temperature dependency. The parameters used in calculating the dielectric function of 


6 



Ini_ x Ga z Sb have been interpolated from InSb and GaSb data using Eq. (I). 5 Cooling down the device 
shifts the dielectric function to higher energy corresponding to shorter wavelength, agreeing with 
reported data. ’ The absorption coefficient, a(A,T), and reflection coefficient, R(A,T), can now be 
calculated using the relations 8 


a(XJ)= T-i 2 


(3) 


r{a,t) 


{a. t ) 2 +e 2 (x, t) 2 - ^2e 1 (A,T)+2j£ l (A,T) 2 +s 2 [ijf + 1 
+ s 2 (A,t) 2 + ^2s 1 (A,T) + 2^ 1 (A,T) 2 +s 2 (Xjf +\ 


(4) 


Fig. 2 shows the absorption and reflection coefficients for Ino. 2 Gao. 8 Sb at temperatures of 253, 273 
and 293 K within the wavelength range of interest to the device analysis. 

Applying eqs (3) and (4), the minority carrier optical generation rate G(x,A,T) along the 

26 

device depth ( x - direction in Fig. 1) is given by applying Beer-Famberf s relation according to" 
A-P 0 JA)-a(A,T)\\-R(A,T)\ 

G(x,A,T ) = — — ", (5) 

h ■ c ■ exp[a( A,T)-x\ 

where P op t(A) is the total optical power incident on the photodetector per unit area. 


2.3 Temperature dependence of the mobilities 

The minority carrier mobilities can be determined by analyzing the scattering mechanisms 
due to the ionized impurities and absorption (or emission) of either acoustical or optical phonons. 
According to Mathiesen’s rule, the mobility of the Ini_ z Ga x Sb alloy is determined from its binary 
compound mobilities using Eq. (1), with inverse mobilities for parameters Q and B and no bowing 
effect. Besides, under low field conditions, majority carrier Hall mobilities and minority carrier 
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drift mobilities are very close to each other 1 1 . Thus, the electron (n) and hole (p) mobilities for both 
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InSb and GaSb were taken from existing Hall mobility data, ’ ' which was fitted to the 
temperature extended Caughey-Thomas empirical model for III-V compounds, described by 9 ’ 1 1 


A 


M( N ’T )<n,p> = P„in,<n,p> + 


max,<n,p> 


f 300 ' 


V 1 J 


7 1 .<»./» 


■p 


min,<n,p> 


( 6 ) 


1 + 


N 


N, 


ref ,<n,p> 


f T \ 


V 300y 


7 2.<n,p> 


where p max and ju min are the carrier mobility values at low and high doping concentrations, 
respectively, N re f is the reference doping concentration for which the mobility reduces to half the 
maximum value, and (j)\s a fitting parameter independent of temperature, for T > 150 K. The positive 
temperature coefficients 6i and d 2 serve to accordingly reduce the lattice limited mobility ju max and 
increase N re f with a rise in temperature. Table II shows the parameters used in the extended 
Caughey-Thomas model to calculate Hall mobilities of electrons and holes for both GaSb and 
InSb. 10,1 1,27 


2.4 Minority carrier lifetimes 

The minority carrier lifetime, r, is comprised of a number of fundamental components owing to 
various recombination mechanisms. Considering the nonradiative Shockley-Read-Hall ( tsrh) and 
Auger ( lAuger) recombination lifetimes, the total minority carrier lifetime is obtained from the inverse 

of the sum of their reciprocals, 31 i.e., z {n p) ~ l = r SRH {n p ^ + r Auger ^ p ~ ] . 

The Shockley-Read-Hall trap-assisted recombination lifetime, caused by imperfections 
within the semiconductor, is expressed as 12,13 


8 



( 7 ) 


1 I m <e,h> 

TsRH ’<i’P> o-jv, -y 3^r 

where cr is the capture cross section of minority carriers and N t is the density of traps, assigned the 
values of 1.5x10' m and 1.17x10“ m' , respectively and kg is the Boltzmann constant. 

Of the possible known Auger recombination mechanisms, the most important for InSb-based 
compounds are between the conduction/valence bands (A-l) and through the conduction/heavy- 
hole/light-hole bands (CHLH or A-7). In the case of In 0 . 2 Ga 0 .gSb, the energy bands of the alloy tend 
towards a GaSb configuration, wherein the magnitude of the spin split-off gap A 0 is significant to the 

Eo energy gap. 1718 Then the Auger recombination through the conduction/heavy-hole/spin split-off 

1 8 

bands (CHSH or A-S) becomes significant. Thus, the total Auger lifetime can be expressed as 


1 


r 


Auger, <n,p> 


n <pO,nO> P<p0,n0> 

2 


1 1 1 

7 — I ] — I / — * 

_P<p0,n0>^Al n <pO,nO> T A7 n <pO,nO> T AS _ 


( 8 ) 


In the above equation, n< p o, n o> and p< p o, n o> are the electron and hole carrier concentrations at 
equilibrium, zai, za 7 , and z A s represent the intrinsic lifetimes for A-l, A-7 and A-S recombination. 
As derived in the work of Beattie and Landsberg 15 , the intrinsic A-l lifetime is given by 


z' AX = 3.8x10' 


• 0 + U m f 2 ■ (1 + 2u m ) • exp f + 2U \} E ° 

V\ l+U m) k B T 


— 1*1*2 | 2 

m„ 


where u m =mjm h and F/F 2 is the overlap integral of the periodic part of the electron wave 

15 16 ... ... 17 is 

function. ’ The intrinsic A-7 lifetime is expressed as ’ 
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2.5 Continuity equations and photocurrent densities 

The distribution of minority carriers in a semiconductor material is governed by the 
continuity equation, which considers the drift, diffusion, photogeneration and recombination effects. 
At steady state under low injection conditions and no external voltage, the one dimensional (along 
the x axis) continuity equations for the n and p regions are: 


D 


n 


<9 2 A n 
dx 2 




( 12 ) 


d +a 

p dx 2 p 



(13) 


where An and Ap are the excess minority carrier densities, D n and D p are the minority carrier 
diffusion coefficients, related to the carrier mobilities by Einstein’s relation. The second order linear 
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inhomogeneous ordinary differential Eqs. (12) and (13) are solved for An and Ap under the following 
boundary conditions 


An(Xj) = 0 , 

(14) 

8An S n 

= " An(0) , 

8x x=0 D n 

(15) 

A p{Xj + w) = 0 , 

(16) 

8*1 

X 

II 

8 

II 

O 

( 17 ) 


In the above S n is the surface recombination velocity of electrons, and w is the depletion region 

26 

depth. The surface recombination velocity is given by 

v the N st , (18) 

where v th e = y]3k B T /m e is the electron thermal velocity and N st is the number of surface trapping 

centers per unit area at the boundary region equal to 2.615x10 m . It should be noted that the 
surface recombination is not considered in eq. (17) due to the use of the ternary substrate, with 
relatively infinite depth (>750 pm). 

26 32 

The depletion width is obtained from” ’ 



where N A and N D are the acceptor and donor concentration and Vbi is the junction built-in potential 

26 32 

given by the equation ’ 



( 20 ) 
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Considering the diffusion lengths L <n p> = J D <n p> r <n p> , the general solutions of Eqs. (12) and (13) 
are given by 


An = A 0 exp 


f \ 
x 


y L nJ 


+ Bg exp 


x 

L 


+ Cg exp -a(A,T )x\. 


n J 


( 21 ) 


A P = Cg\ 
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:p -a(A,T )x\ - exp 
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(Xj+W) 
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a(A,T) 

\ L p 


J 


x 

L„ 


( 22 ) 


where 


„ P v (A)\\-R(X,T)\A-a(X,T) 

^ A / \ 


Dhc\ 


a(A,T ) 2 - 


1 


L 


n J 


(23) 


Bg =-^[l + L n a( A,T )\ exp 

n 

A„= ~ ‘ [B 0 (S„L,+D,) + C a L,(D,aa,T) + Sj. 


A-a(Z.T) 

\ L n 


(24) 

(25) 


Once the excess minority carriers distribution are known, the diffusion current densities in 
the p and n regions can be obtained using the relations 


J n = C 1 D , 


GAn 

dx 


(26) 


J n = ~qD_ 


GAp 


dx 


(27) 


for the electron and holes respectively, while the photocurrent density in the depletion region is 
given by 


X j+W 

J scr =-q J G(. x)dx , 


(28) 
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The total photocurrent density for the device J to t is the sum of the individual contributions from each 
region, i.e., J tot = J n + J scr + J p . By definition, the responsivity, Di, is then calculated from 


M(A,T) 


p oP ,(v 


and therefore the quantum efficiency, rj, is obtained from 

J_JM) he 

q K r ,W 
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(29) 


(30) 


3. RESULTS AND DISCUSSION 

An Ino. 2 Gao. 8 Sb p-n junction photodetector was formed by Zn-diffusion to an Ino. 2 Gao. 8 Sb 
virtual substrate using the leaky box technique 20 . The initial substrate doping was Ad=1x10 17 cm" 3 . 
The diffusion depth was Xj=0A pm, with acceptor concentration of Na= 1x10 cm" . The 
characterization results of this device were presented elsewhere" . Fig. 3 shows the spectral response 
of a 300x300 pm 2 InGaSb photodetector compared to the model results at three different 
temperatures. The corresponding quantum efficiency as obtained from Eq. (4) is presented in Fig. 4. 
The cut-off wavelength was observed to shift to a shorter value with an increase in the responsivity 
and quantum efficiency when cooling the device. The shift in the cut-off wavelength is directly 
related to the temperature dependence of the dielectric function, through the energy band gaps, 
presented in both absorption and reflection coefficients. The increase in the responsivity and 
quantum efficiency at lower temperatures is due to the variation of minority carrier recombination. 
Reducing the device temperature leads to suppress thermal vibrations, resulting in increase in the 
minority carriers mobility and lifetime. Thus, the diffusion length increases leading to the observed 
increase in the perfonnance. According to the model, for the electrons in the diffused p-region, SRH 
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is the dominant recombination mechanism, which has relatively lower temperature dependence. For 
the n-region, presented by the ternary substrate, Auger recombination dominates, especially A-S 
with relatively higher temperature dependence, down to approximately 180 K. Below that 
temperature, SRH becomes dominant. 

Focusing on the room temperature spectral response (293 K), the inset of Fig. 4 shows the 
individual contribution from each region to the total responsivity of the device. As indicated in the 
figure the n region has higher influence on the performance compared to the regions at longer 
wavelengths. This is an advantage due to the ternary substrate compared to a binary one. The 
responsivity decreases at shorter wavelengths due to the domination of surface recombination. When 
the photon energy is less than the energy band gap, the responsivity tends to zero. For this device, 
this occurs at approximately 2100 nm. At lower photon energies, most carriers are generated in the n 
region since the absorption length is longer. As the photon energy increases, the absorption length is 
reduced and the responsivity of the p region becomes dominant. This occurs at a wavelength of 1 100 
nm for this device. Focusing on this fact, Fig. 5 presents the effect of the Zn-diffusion depth, x h on 
the quantum efficiency. At shorter wavelength, increasing the diffusion depth (the p-region depth) 
increases the recombination probability of photogenerated carriers, thus reducing the quantum 
efficiency. 

The effect of Zn-diffusion depth and temperature on responsivity is future discussed in the 
inset of Fig. 5, assuming an arbitrary wavelength of 1.9 pm. As the diffusion depth increases, fewer 
photogenerated carriers from the p-region reach the depletion layer, thereby reducing the 
responsivity. At a given temperature, if the length is shorter than a certain value (optimum diffusion 
depth), the responsivity decreases because fewer photons can get absorbed in such a narrow p- 
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region. The optimum diffusion depth has been found to decrease with cooling of the device, with a 
value of 0.3 pm at 293 K. 

The responsivity for various acceptor and donor concentrations for the same Ino. 2 Gao. 8 Sb 
photodetector is shown in Fig. 6. Increasing the concentration of acceptors and/or donors slightly 
augments the built-in potential that serves to separate the generated electron-hole pairs (Eq. 20). 
However, the overall effect is a decrease in the depletion region width, where the photogeneration of 
carriers mainly takes place (Eq. 19). Thus, higher doping levels adversely affect responsivity, as can 
be seen in the figure. On the other hand, very low doping levels are not practically achievable 
besides it deteriorates the transient response. 

4. CONCLUSIONS 

The spectral response of In o. 2 Ga o.sSb p-n junction infrared photododetectors has been 
analyzed using a theoretical model based on the continuity equations of the minority carriers in the 
p, n and space charge regions. The absorption coefficient and normal-incidence reflectivity have 
been incorporated using a dielectric function model developed by Adachi and extended here to 
account for temperature effects. The carrier mobility has been calculated using the temperature- 
extended Caughey-Thomas model. The changes in the absorption coefficient and reflectivity with 
decreasing temperature cause a shift in the cutoff wavelength to shorter values while the Auger 
recombination mechanism have proven to have a dominating effect in the horizontal shift of the 
spectral response with cooling. These effects correspond to the energy bands variation with 
temperature, which has been approximated using Varshini’s equation. The effects on the spectral 
response due to the variation of temperature, photon energy, diffusion depth, and carrier 
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concentration have all been analyzed. By appropriate combination of the aforementioned parameters, 
the performance of the device can be optimized. 
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Table I. Temperature dependence of the critical-point energies of GaSb and InSb, with parameters 
used in Varshni’s equation for each energy level (after Refs. 23-25). 


Critical Point 
Energies 


GaSb 



InSb 


E(0) 

(eV) 

8 

(xlO' 4 eV.K' 1 ) 

P 

(K) 

E(0) 

(eV) 

8 

(xlO' 4 eV.K' 1 ) 

P 

(K) 

E 0 

0.812 

4.2 

140 

0.24 

6 

500 

E 0 +A 0 

1.6 

5 

143 

1.2 

6 

500 

E, 

2.186 

6.8 

147 

2 

6.84 

132 

Ei+A] 

2.621 

6.7 

176 

2.49 

6.46 

170 

e 2 

4.32 

9 

376 

4.24 

5.4 

0 
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Table II. Parameters used in the temperature extended Caughey-Thomas model for calculating the 
electron (e) and hole (h) mobilities for GaSb (Ref. 11) and InSb (Ref. 27). 


Parameters 


Material/Carrier 


GaSb/e 

GaSb/h 

InSb/e 

InSb/h 

Umax (cm 2 /V-s) 

5650 

875 

7.8xl0 4 

750 

JUmin (Cm 2 /V-S) 

1050 

190 

5000 

100 

yV re /(cm" 3 ) 

2.8xl0 17 

9xl0 17 

7xl0 16 

6xl0 17 

4> 

1.05 

0.65 

0.7 

0.6 

9i 

2.0 

1.7 



e 2 

2.8 

2.7 
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List of Figure Captions 


FIG. 1 Schematic cross section of the modeled Ini_ x Ga x Sb p-n photodetector, xj is the Zn-diffusion 
depth and w is the depletion region width. 

FIG. 2 Ino. 2 Gao. 8 Sb absorption coefficient and normal-incidence reflection coefficient versus 
wavelength calculated at temperatures of 253, 273 and 293 K. 

FIG. 3 Spectral response at different temperatures obtained from the model compared to the 
experimental data. The model results were obtained assuming diffusion depth xj = 0.4 pm, and 
acceptor and donor concentrations of N A = lxlO 19 cm" 3 and N D = lxlO 17 cm" 3 , respectively. 

FIG. 4 Calculated quantum efficiency variation with wavelength at different temperatures. The 
results obtained using the same parameters of Fig. 3. The inset shows the 293 K calculated 
responsivity of the different regions of the device compared to the total responsivity. 

FIG. 5 The calculated variation of the quantum efficiency with the Zn-diffusion depth, obtained at 
293 K. Focusing on 1.9 pm wavelength, the inset shows the responsivity variation with the Zn- 
diffusion depth at three different temperatures. 

FIG. 6 Responsivity variation with Zn-diffusion concentration (acceptor level Na) and the original 
doping of the ternary substrate (donor level No). Calculations assume 0.8 composition, 293 K 
temperature, 0.4 pm Zn-diffusion depth and 1.9 pm wavelength of the incident radiation. 
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